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We report the measurement of near threshold neutral pion electroproduction cross sections and 
the extraction of the associated structure functions on the proton in the kinematic range Q 2 from 
2 to 4.5 GeV 2 and W from 1.08 to 1.16 GeV. These measurements allow us to access the dominant 
pion-nucleon s-wave multipoles Eo+ and So+ in the near-threshold region. In the light-cone sum- 
rule framework (LCSR), these multipoles are related to the generalized form factors G\ P (Q 2 ) and 

P {Q 2 )- The data are compared to these generalized form factors and the results for GJ P (Q 2 ) 
are found to be in good agreement with the LCSR predictions, but the level of agreement with 
Gj P (Q 2 ) is poor. 

PACS numbers: 25.30.Rw, 13.40.Gp 



I. INTRODUCTION 

Pion photo- and electroproduction on the nucleon 
jN — > ttN, j*N —> ttN close to threshold has been 
studied extensively since the 1950s both experimentally 
and theoretically. Exact predictions for the threshold 
cross sections and the axial form factor were pioneered 
by Kroll and Rudcrman in 1954 for photo-production and 
are known as the low energy theorem (LET) jT] . This 
LET provided model independent predictions of cross 
sections for pion photoproduction in the threshold re- 
gion by applying gauge and Lorentz invariance [2 . This 
was the first of the LET predictions to appear but was 
not without limitations. This LET predictions were re- 
stricted only to charged pions and the it contribution 
was shown to vanish in the 'soft pion' limit, i.e., ~ p v . 
Here, andp^ are the mass and momentum of the pion. 
Additionally, these cross section predictions were limited 
to diagrams with first order contributions in the pion- 
nucleon mass ratio. In later years, using vanishing pion 
mass chiral symmetry (m n — > 0), these predictions were 
extended to pion electroproduction for both charged and 
neutral pions [21 H] • 

Of course, a vanishing pion mass doesn't relate to the 
observed mass of the pion (the pion to nucleon mass ratio 
m^/rriM ~ 1/7), so higher order finite mass corrections 
to the LET were formulated in the late sixties and early 
seventies before the appearance of QCD. These also in- 
cluded contributions to the non-vanishing neutral pion 
amplitudes for the cross section. 

In the late eighties and early nineties, experiments at 
Mainz 5j obtained threshold pion photo-production data 



on 7p — > ir°p. The theoretical predictions of LETs at 
the time were inconsistent with the data at low pho- 
ton energies. With the emergence of chiral perturba- 
tion theory (%PT), the scattering amplitudes and some 
physical observables were systematically expanded in the 
low energy limit in powers of pion mass and momen- 
tum. Using this framework, the LET was re-derived to 
include contributions to the amplitudes from certain loop 
diagrams, which were lost when the expansion was per- 
formed in terms of the pion mass, as was done in the ear- 
lier works [7] . Further electroproduction experiments 
at NIKHEF [8] on 7*p -> ir°p with photon virtuality 
Q 2 ~ 0.05 — 0.1 GeV 2 [9] provided good agreement with 
%PT predictions. 

These LETs [B H II El H] are not applicable for 
Q 2 > A 3 QCD / m7r , where A QC d ~ 200 - 300 MeV is 
the QCD scale parameter. In the case of asymptotically 
large momentum transfers (Q 2 — > oo) perturbativc QCD 
(pQCD) factorization techniques [TUMIi?] have been used 
to obtain predictions for cross section amplitudes and 
axial form factors near threshold. In these factorization 
techniques, 'hard' (Q 2 3> Aq CD ) and 'soft' (fc ~ Aqcd) 
momentum contributions to the scattering amplitude can 
be separated cleanly and each contribution can be theo- 
retically calculated using pQCD and LETs, respectively. 
Here, k is the momentum of the virtual photon. 

Recently, Braun et al. [121 [T3] suggested a method 
to extract the generalized form factors, G\ N (Q 2 ) and 
Gl N {Q 2 ), for 1 < Q 2 < 10 GeV 2 using light cone sum 
rules (LCSR). The transition matrix elements of the elec- 
tromagnetic interaction, J M , can be written in terms of 
these form factors at threshold: 
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<iV(P>(fc)| J>(P)) = - t n 15 
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(Q 2 ) 
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Here, N(P') and p(-P) are spinors for the final and initial 
nucleons with momenta P' and P, respectively, mjv is 
the mass of the nucleon, f v is the pion decay constant 
and q is the 4-momentum of the virtual photon. Since the 
pion is a negative parity particle and the electromagnetic 
current is parity conserving, the 75 matrix is present to 
conserve the overall parity of the reaction. 

These form factors are directly related to the pion- 
nucleon s-wave multipoles Eq + and L$ + 03 HI] 
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Here, a is the electromagnetic coupling constant and 
ui!y h is the virtual photon energy at threshold in the 
cm. frame and is given by the following relation: 

2(771 JV + "V) 



th _ m w (2m N + m w ) - 

W 7 - 



(4) 



In general, E\±, Mi±, and L;± describe the electric, mag- 
netic and longitudinal multipoles, respectively. Here, I 
describes the total orbital angular momentum of the pion 
relative to the nucleon and ± is short for ±| so that the 
total angular momentum of the ttN system is I ± , • 

Additionally, the sum rules can be extended to the 
Q 2 ~ 1 GeV 2 regime and the LETs are recovered to 
0(m„-) accuracy by including contributions from semi 
disconnected pion- nucleon diagrams [14] . This approach 
provides a connection between the low and high Q 2 
regimes. Predictions for the axial form factor and the 
generalized form factors are also obtained in this ap- 
proach. 

In the low Q 2 < 1 GeV 2 regime and the chiral limit 
TOtt — > 0, the LET s-wave multipoles at threshold can be 
written as 171: 



Eq+ — 



8tt 



G\ 
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(5) 
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Gi and GJ can be written in terms of the electro- 
magnetic form factors for the neutral pion-proton 7r°p 
channel in this approximation: 



Q 2 



G\ 



fjA 
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(7) 
(8) 



In the above equations, G P M and G P E are the Sachs elec- 
tromagnetic form factors of the proton and gA is the axial 
coupling constant obtained from weak interactions. Also, 
for the charged pion- neutron 7r + n channel, the general- 
ized form factors can be written as: 
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2m 2 N 



G\ 



(10) 



Here, G\ t and G% are the electromagnetic form factors 
of the neutron. Additionally, G a is the axial form factor 
that is induced by the charged current and its contribu- 
tion comes from the Kroll-Ruderman term pQ. 

These generalized form factors, GJ N and GJ W , can 
be described as overlap integrals of the nucleon and the 
pion-nucleon wave functions. The wave function of the 
pion-nucleon system at threshold is related to the nu- 
cleon wave function without the pion by a chiral rotation 
in the spin-isospin space [lOj [13] . The measurement of 
these form factors for pion clectroproduction is in essence 
the measurement of the overlap integrals of the rotated 
and non-rotated nucleon wave functions, which are not 
accessible in elastic form factor measurements. This in- 
formation complements our understanding of the various 
components of the nucleon wave function (quarks and 
gluons) and the theory of strong interactions. Addition- 
ally, it provides insight into chiral symmetry and its vio- 
lation in reactions at increasing Q 2 . 

The generalized form factor for the charged pion- 
neutron G\ n (Q 2 ) and the axial form factor Ga{Q 2 ) had 
been measured near threshold for Q 2 ~ 2 — 4.2 GeV 2 [15j . 
In this paper, we describe the measurement of the dif- 
ferential cross sections and the extraction of the s-wave 
amplitudes for the neutral pion electroproduction pro- 
cess, ep —¥ epir , for Q 2 ~ 2 — 4.5 GeV 2 near threshold, 
i.e., W ~ 1.08— 1.16 GeV. From these cross sections, the 

generalized form factors G\ P (Q 2 ) and G\ P {Q 2 ) were ex- 
tracted and compared with the theoretical calculations of 
Refs. HH and [7J. 
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FIG. 1. Neutral pion electroproduction in the center of mass 
frame. 



II. KINEMATIC DEFINITIONS AND 
NOTATIONS 

The neutral pion reaction 

e(0 + P (P) -> e(0 + p{P') + 7r°(fc) (11) 

is shown schematically in the virtual photon-proton cen- 
ter of mass frame in Fig. [I] Here, I = (E e ,p e ), V — 
(E' e ,p' e ), P = (m p ,0) and P' = (E' p ,p' p ) are the ini- 
tial and final electron and proton 4-momenta in the lab 
frame and k = (E„, p w ) is the 4-momentum of the emit- 
ted pion. Also, m p refers to the mass of the proton. It is 
assumed that the incident electron interacts with the tar- 
get proton via exchange of a single virtual photon with 
4-momentum q = I — I' = (u), q). In this approximation, 
it is also assumed that the electron mass is negligible 
(m e ~ 0). The two important kinematic invariants of 
interest are 



Q 



4E e E' e sm 2 (9' e /2) 



s = W 2 = {q + Pf 



m p + 2u)m v 



Q ■ (12) 



Here, 9' e is the polar angle of the scattered electron in the 
lab frame. 

The five-fold differential cross section for the reaction 
can be written in terms of the cross section for the sub- 
process 7*p — > pir° [TS] , which depends only on the ma- 
trix elements of the hadronic interaction: 



d 5 a 



dE'dSKdSli 



cr -/*p 
dSil 



(13) 



Here, dd' e — d cos B' e d<j>' e is the differential solid angle 
for the scattered electron in the lab frame and dCl* — 
dcos6*d<f>* is the differential solid angle for the pion in 
the virtual photon-proton {"f*p) center of mass frame. 
The azimuthal angle </>* is determined with respect to 
the plane defined by the incident and scattered lepton 
[2J. The factor T represents the virtual photon flux. In 
the Hand convention fTU it is 



a 



E' W 2 



1 



2tt 2 E e 2m p Q 2 1-e' 



(14) 



which depends entirely on the matrix elements of the lep- 
tonic interaction and contains the transverse polarization 
of the virtual photon 



2 |q|2 tan 2 ^ 



(15) 



For unpolarized beam and target the reduced cross 
section from Eq. ( 13 ) can be expanded in terms of the 



hadronic structure functions: 
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(16) 



Here, p* is the pion momentum and K = (W 2 — m p ) /2W 
is the photon equivalent energy in the cm. frame of the 
subprocess j*p —> pir°. Additionally, ot cflt and 

ott are the structure functions that describe the trans- 
verse, longitudinal, longitudinal-transverse interference, 
and transverse-transverse interference components of the 
differential cross section. 

Each of these structure functions contain the cos6>* 
dependence and can be parameterized in terms of the 
multipole amplitudes Ei±, Mj± and Si± that describe 
the electric, magnetic and scalar multipoles, respectively. 
The scalar multipoles Si± can be written in terms of the 
longitudinal multipoles Li± — j^Si±, where lu* and q* 
are the energy and 3-momentum of the virtual photon in 
the cm. frame, respectively [2]. 



III. EXPERIMENT 

The near threshold reaction ep — > epir was studied us- 
ing the CEBAF Large Acceptance S pectro meter (CLAS) 
in Jefferson Lab's Hall-B [T7]. Fig. 2(a)| shows the de- 
tector components that comprise CLAS. Six supercon- 
ducting coils of the torus divide CLAS into six identical 
sectors and produce a toroidal magnetic field in the az- 
imuthal direction around the beam axis. Each of the six 
sectors contain three regions of drift chambers (Rl, R2, 
and R3) to track charged particles and to reconstruct 
their momentum [18] . scintillator counters for identify- 
ing particles based on time-of-flight (TOF) information 
|19) . Cerenkov counters (CC) to identify electrons (20) . 
and electromagnetic counters (EC) to identify electrons 
and neutral particles [3T]. The CC and EC are used for 
triggering on electrons and provide a mechanism to sepa- 
rate charged pions and electrons. With these six sectors, 
CLAS provides a large solid angle coverage with typical 
momentum resolutions of about 0.5% — 1.0% depending 
on the kinematics [T7] . 

A 5.754 GeV electron beam with an average intensity 
of 7 nA was incident on a 5 cm long liquid hydrogen tar- 
get, which was placed 4 cm upstream of the CLAS center. 
Fig. 2(a) shows the electron beam entering CLAS from 
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(a) 



3 3.5 4 4.5 
p (GeV) 

(a) 




FIG. 2. (a) A three-dimensional view of CLAS showing the 



superconducting coils of the torus, the three regions of drift 
chambers (R1-R3), the Cerenkov counters, the time-of- flight 
system, and the electromagnetic calorimeters. The positive 
z-axis is out of the page along the symmetry axis, (b) A 



schematic view of a typical near threshold event showing the 
reconstructed electron and proton tracks with the correspond- 
ing detector hits in two opposite CLAS sectors. The 7T° is 
reconstructed using the missing mass technique as discussed 
in the text. 



the top left and exiting from the bottom right through 
the symmetry axis. A small non-superconducting mag- 
net (minitorus) surrounded the target and generated a 
toroidal field to shield the Rl drift chambers from low 
energy electrons of high intensity. These electrons origi- 
nated primarily from the M0ller scattering process. The 
data used in this experiment were collected from Octo- 
ber 2001 to January 2002 and the integrated luminosity 
was about 0.28 fb -1 . The electron beam energy of 5.754 
GeV as determined in this experiment agrees within 6 
MeV with an independent measurement in Hall A [22] . 



IV. ANALYSIS 

At the start of this analysis, a cut of W < 1.3 GeV 
is applied to focus our events only in the kinematic re- 




(b) 



FIG. 3. (Color online) EC sampling fraction as a function of 
electron momentum for one of the CLAS sectors for (a) Data 
and |(b)| Monte Carlo (MC) simulation. The dashed lines show 
the parameterized mean and the solid line indicates the 3a 
cut. 



gion of interest. In this analysis the scattered electrons 
and protons are detected using CLAS and the ir° is re- 
constructed using 4-momentum conservation. A typical 
event for this experiment is shown in Fig. |2(b)| 



A. Particle Identification: Electron 

The scattered electrons in the final state of the reac- 
tion are detected by requiring geometrical coincidence 
between the Cerenkov counters and the electromagnetic 
calorimeter in the same sector. The momentum of the 
electrons is reconstructed using the drift chambers. Using 
the energy deposited in the EC and the momentum, the 
electrons are isolated from most of the minimum ionizing 
particles (MIPs), e.g., pions, contaminating the electron 
spectra. 

As electrons pass through the EC, they shower with a 
total energy deposition E to t that is proportional to their 
momenta p. The sampling fraction energy E tot /p is plot- 
ted as a function of momentum for each sector after ap- 
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i.e., 




(a) 




(b) 



FIG. 4. (Color online) At as a function of p. The curves show 
the ±3.5cr cut (solid lines) from the mean fit (dashed line) for 
one of the CLAS sectors for |(a)| experimental and |(b)| Monte 
Carlo simulated events. 



plying all the other electron identification cuts. Fig. [3] 
shows this distribution for one of the CLAS sectors for 
experimental and Monte Carlo simulated events. In the 
figure, one can note the MIPs contamination near the 
smaller values of E to t/p- This contamination is signifi- 
cantly larger in data than in simulated events. The elec- 
trons are concentrated near E to t/p ~ 0.3. Ideally they 
should not show any dependence on momentum, albeit a 
slight momentum dependence is visible in the data. This 
dependence is parameterized and a cut of 3cr is applied as 
shown in the figure. The MIP events are well separated 
from the electrons below the 3a cut. 



B. Particle Identification: Proton 

The recoiled protons are identified using the measured 
momentum and the timing information obtained from 
the TOF counters. A track is selected as a proton whose 
measured time is closest to that expected of a real proton, 



At = t r , 



(t 



TOF 



0c 



dcC 



(17) 



In the above equation, txoF is the time measured from 
the TOF counters, I is the distance from the target cen- 
ter to the TOF paddle, and t tr is the event start time 
calculated from the electron hit time from the TOF 
traced back to the target position. Also, in Eq. (17) 

Pcaic = Pl\J M ldg where /3 ca/c is computed using 

the PDG [23] value of the mass of the proton M p d g and 
the momentum of the track p. 



Figs. 4(a) and |(b)| show the experimental and simulated 
event distributions, respectively, of At as a function of p 
for one of the CLAS sectors. The protons are centered 
around At = ns and have a slight momentum depen- 
dence for p < 1 GeV. The dashed lines indicate the pa- 
rameterized mean of the distributions and the solid lines 
indicate the ±3.5ct cut applied to select the protons. 



C. Fiducial Cuts and Kinematic Corrections 

For perfect beam alignment, the incident electron 
beam is expected to be centered at (Xb ea mi Ybeam) — 
(0, 0) cm at the target. But due to misalignments, 
the electron beam was actually at (Xb eam ,Yi, eam ) = 
(0.090, —0.345) cm. This misalignment of the beam- 
axis is corrected for each sector, which also subsequently 
changes the reconstructed z- vertex positions of the elec- 
tron and proton tracks. The details of this correc- 
tion are described in previous works [24] [25] . A cut of 
z € (—8.0, —0.8) cm is placed on the z-vertex to isolate 
events from within the target cell. 

The measured angles and momenta of the electrons and 
protons are corrected using the same method as used in 
previous analyses [2"4"1 12"5] . 

The electrons start to lose energy as they enter the 
electromagnetic calorimeter. When the electrons shower 
near the edge of the calorimeter, their shower is not fully 
contained and so their energies cannot be properly re- 
constructed. As such, a fiducial cut is applied to remove 
these events. 

Electrons give off Cerenkov light in the CC, which is 
collected in the PMTs on either side of the counters in 
each sector. Inefficient regions in the CC are isolated 
by removing those regions where the average number of 
photo-electrons (Nphe) < 5. This cut results in keeping 
all events that lie in regions where the CC efficiency is 
about 99% [2D]. 

To deal with edges and holes in the drift chambers, and 
to remove dead or inefficient wires, a fiducial cut for both 
electrons and protons is applied. Regions of non-uniform 
acceptance in the azimuthal angle <j> resulting from these 
attributes are isolated on a sector-by-sector basis as a 
function of the electron's momentum p e and polar angle 
e . For the electron, at fixed p e and 9 e , one expects the 
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FIG. 5. Electron cf> e distribution for CLAS Sector 4 for 
p e — 4.1 ±0.1 GeV shown for different e slices. The unshaded 
curves show <f) e distribution after electron selection and the 
shaded curves show the <f) e distribution after applying electron 
DC fiducial cuts. 
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FIG. 6. (Color online) Proton cj> p vs. 6 V distribution for CLAS 
Sector 4 for p p = 2.85 ± 0.15 GeV. Rejected tracks are shown 
in black. 



angular distribution to be symmetric in <f> e and relatively 
flat. Empirical cuts are applied to select these regions 
of relatively flat <f> e as shown in Fig. [5] for electrons with 
p = 4.1 ± 0.1 GeV for different slices of 8 e and one of 
the CLAS sectors. The same cuts are applied to both 
experimental and simulated events. 

As for electrons, a fiducial cut on the proton's az- 
imuthal angle <j> p as a function of its momentum p p and 
polar angle 8 p is applied. However, the edges of the <j)p 
distributions are asymmetric for different slices of 6 p . 
The upper and lower bounds on <j> p are extracted and 
parameterized as a function of 9 p and p p . The result of 
this cut for one of the CLAS sectors is shown in Fig. [6| 
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FIG. 7. The Bethe-Heitler process ep — > ep-y diagrams for (a) 
a photon emitted from an incident electron (pre-radiation) 
and for (b) a photon emitted from a scattered electron (post- 
radiation) - ^ 



D. Background Subtraction and ir° Identification 

The neutral pion in the final state is reconstructed us- 
ing energy and momentum conservation constraint. To 
do so, we use the conservation of 4-momentum and look 
at the missing mass squared distribution of the detected 
particles (i.e., the electron and the proton): 



M 2 x {ep) = {l 



P-V- P') 



(18) 



Here, I, P, I' and P' are 4-momenta of the incident and 
scattered particles as described in Section [Tl] 

There are several difficulties in the analysis in the near 
threshold region. In this region, the pion electroproduc- 
tion cross section goes to zero; so, the statistics are very 
low. Also, a major source of contamination to the neutral 
pion signal near threshold is the elastic Bethe-Heitler pro- 
cess ep — > epj. The two dominating Feynman diagrams 
for this process are shown in Fig. [7j Fig. 7(a) shows 



the diagram with a pre-radiated photon (emission from 
an incident electron) and Fig. |7(b)| shows the diagram 
with a post-radiated photon (emission from a scattered 
electron). These photons are emitted approximately in 
the direction of the incident and scattered electron, re- 
spectively H7] . When these photons are emitted, the 
incident and scattered electrons lose energy. This fea- 
ture of the Bethe-Heitler process can be exploited to our 
benefit. 

For the elastic process ep — > ep, the proton angle can 
be computed independently of the incident or scattered 
electron energies: 



tan0f 



tan #2 = 



1 



E 1 



rrip — E' cos t 
1 



tan -£ 



(19) 
(20) 



tan -f 



Here, Q\ and Q\ are the proton angles computed inde- 
pendently of the incident or scattered electron energies, 
respectively. Also, Q' e is the angle of the scattered elec- 
tron in the lab frame, and E and E are the energies 
of the incident and scattered electron, respectively. We 
can calculate these angles for each event and look at its 
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FIG. 8. (Color online) | (a) | Af| vs A6{ for W = 1.09 ± 0.01 
GeV. The red dashed line indicates the expected pion peak 
position. The left red spot centered around zero degrees cor- 
responds to the elastic scattering events in which the inci- 
dent electrons have undergone Bethe-Heitler radiation (pre- 
radiative) and the one on the right to the elastic post-radiative 
events. The events below the linea r po lynomial and outside 

My for events with 



(b) 



the ellipse are selected as pions. 
W = 1.09±0.01 GeV. The black solid curve shows events prior 
to any Bethe-Heitler subtraction cuts, the blue dashed-dot 
curve shows events rejected from the cuts, and the red dashed 
curve shows those events that survive the Bethe-Heitler sub- 
traction cuts. 



deviation (A6^ 2 ) from the measured value {O^eas)'- 



A0? >2 = 91 



1,2 



(21) 



Fig. 



shows the M\ plotted as a function of this 
~A6\ for one of the near threshold regions, 
GeV. In the plot, we see two red 
: OGeV 2 . The one on the left is 
deg corresponding to the pre- 
The other corresponds to the 
Additionally, these radia- 
Thcse are 



deviation 

W = 1.09 ± 0.01 

spots along M\ — 
centered along AQ\ 
radiated photon events, 
post-radiated photon events, 
tive events are also present in the positive M\ . 
the radiative events that we need to isolate from the pion 
signal as indicated by the red dashed line in the plot. An 
ellipse and a linear polynomial are used to reject these 
events. These cuts are parameterized as a function of W. 
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FIG. 9. (Color online) An example of the M%(ep) distri- 
bution with a double Gaussian fit after applying the ellipti- 
cal cuts (black circles) of Fig. 8(a) and after residual Bethe- 



Heitler and other contamination subtractions (green trian- 
gles) for Q 2 = 2.75 ± 0.25 GeV 2 and W = 1.09 ± 0.01 GeV 
(top) and W = 1.11 ± 0.01 GeV (bottom) integrated over all 
<j>Z and cos#£. The black dashed lines indicate the ±3<r cuts 
applied to select the pions. The \ 2 i s the goodness of fit per 
degree of freedom. See Sec. |IV D] for details. 



The result of these cuts is seen in Fig. 8(b) with the ac- 
cepted events after the cut shown in red (dashed curve) 
as our pions and the rejected events in blue (dashed-dot 
curve) . 

After the Bethe-Heitler subtraction cuts are applied, 
the pions are selected by making a ±3ct cut on M\ from 
the mean position of the distribution. An example of the 
distributions and fit are shown in Fig. [9j The M\ distri- 
butions (black circles) are fit with two Gaussians. The 
blue (dashed-dot) curve is an estimate of the remaining 
Bethe-Heitler background in the M\ distribution, which 
was not eliminated by the elliptical cuts of Fig 



This was subtracted to yield the green (triangle) points. 
A systematic uncertainty of ±8% is associated with this 
background subtraction procedure, which is detailed in 
Sec. En 



V. SIMULATIONS 

To determine the cross section, a Monte Carlo simu- 
lation study is required, including a physics event gen- 
erator and the detector geometry. Events are generated 
using the MAID2007 unitary isobar model (UIM) [3H], 
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Variable 


Range 


Number of Bins 


Width 


W (GeV) 


1.08 : 1.16 


4 


0.02 


Q 2 (GeV 2 ) 


2.0 : 4.5 


4 


variable 


cos #* 


-1 : 1 


5 


0.4 


<t>l (deg) 


: 360 


6 


60 



TABLE I. Kinematic bin selection. 



which uses a phenomenological fit to previous photo- 
and electroproduction data. Nucleon resonances are de- 
scribed using Breit-Wigner forms and the non-resonant 
backgrounds are modeled from Born terms and t-channel 
vector-meson exchange. To describe the threshold behav- 
ior, Born terms were included with mixed pseudovector- 
pseudoscalar n N N coupling [5S] . While the pion electro- 
production world-data in the resonance region goes up to 
Q 2 ~ 7 GeV 2 [29] for W > 1.11 GeV, there are no data 
near threshold for Q 2 > 2 GeV 2 and W < 1.11 GeV (the 
kinematics of this work). Thus, cross sections for the 
kinematics of this work are described by extrapolations 
of the fits to the existing data in the MAID2007 model. 

Events are generated to cover the entire kinematic 
range described in Table |TJ About 73 million events are 
generated for the 2400 kinematic bins and 6.7 million 
events were reconstructed after all analysis cuts. The 
average resolutions of the kinematic quantities, W, Q 2 , 
cos6*, and 0; are 0.014 GeV, 0.008 GeV 2 , 0.05, and 8 
degrees, respectively. These resolutions are obtained by 
comparing the generated kinematic quantities with those 
after reconstruction. 

After the physics events are generated, their passage 
through the detector is simulated using the GEANT3 
based Monte Carlo (GSIM) program. This program sim- 
ulates the geometry of the CLAS detector during the 
experiment and the interaction of the particles with the 
detector material. GSIM models the effects of multiple 
scattering of particles in the CLAS detector and geomet- 
ric mis-alignments. The information for all interactions 
with the detectors is recorded in raw banks, which is used 
for reconstruction of the tracks. 

The events from GSIM are fed through a program 
called the GSIM Post Processor (GPP) to incorporate 
effects of tracking resolution and dead wires in the drift 
chambers, and timing resolutions of the TOF. 
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FIG. 10. Acceptance corrections for W = 1.09 GeV and 
Q 2 = 2.75 GeV 2 as a function of </>*.. Each subplot shows the 
correction for a different cos (9* bin. 



VI. CORRECTIONS 
A. Acceptance Corrections 

Acceptance corrections are applied to the experimen- 
tal data to obtain the cross section for each kinematic 
bin. These corrections describe the geometrical coverage 
of the CLAS detector, inefficiencies in hardware and soft- 
ware, and resolution effects from track reconstruction. 

By comparing the number of events in each kinematic 
bin from the physics generator and the reconstruction 
process, the acceptance can be obtained as: 



A 



NL 



N, 



(22) 



yen 



where N^ ec corresponds to those events that have gone 
through the entire analysis process including track re- 
construction and all analysis cuts. N l gen are those events 
that were generated. Fig. [10| shows the acceptances for a 
few of the near threshold bins as a function of 6* . 



B. Radiative Corrections 

The radiative correction is obtained using the software 
package EXCLURAD [3U] that takes theoretical models 
as input to compute the corrections. For this experiment 
the MAID2007 model, the same model used to generate 
Monte Carlo events, is used to determine the radiative 
corrections. The radiative corrections are closely related 
to the acceptance corrections. For each kinematic bin the 
differential cross section can be written as: 



These events are then processed using the same codes 
as those events from the experiment to reconstruct tracks 
and higher level information such as 4-momentum, tim- 
ing, and so on. The simulated events are analyzed the 
same way as the experimental data and are used to ob- 
tain acceptance corrections and radiative corrections for 
the cross sections calculations. 



LA <T 



(23) 



where N meas /C is the number of events from the ex- 
periment normalized by the integrated luminosity (with 
appropriate factors) before acceptance and radiative cor- 
rections. Also, A = N^ D /N^ D is the acceptance cor- 
rection for the bin and 6 is the radiative correction. It 
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W: 1.11 GeV 

: 3.25 GeV 2 




FIG. 11. The radiative corrections for W = 1.11 GeV and 
Q 2 = 3.25 GeV 2 as a function of cos 9* and obtained from 
EXCLURAD using the MAID2007 model. 



should be noted that the events for the acceptance cor- 
rection were generated with a radiated photon in the final 
state using the MAID2007 model. 

EXCLURAD uses the same model to obtain the correc- 
tion S = N£i D ' /N£° RAD ' , where N™ RAD ' are events 
generated without a radiated photon in the final state. 
Thus 



/ atRAD \ ( atNORAD' ' 
meas I gen \ I gen 

~T I NRAD I X I N RAD' 



(24) 



The details of the radiative correction procedure are de- 
scribed in Ref. [2"5] . 

Fig. [TT] shows the radiative corrections calculated for 
one of the kinematic bins as a function of the pion an- 
gles in the cm. system. One can observe that the cor- 
rections have a dependence. This is because the 
bremsstrahlung process only occurs near the leptonic 
plane, i.e., at angles near or 180 degrees with respect 
to the hadronic plane. Also, one can notice that the cor- 
rection increases with cos 0* — > — 1. This is because the 
cross section is expected to approach zero at backwards 
angles and that is the region where the Bethe-Heitler 
events dominate. The average radiative correction over 
all kinematic bins is ~ 25%. 



C. Other Corrections 

Two other corrections were applied to the cross section. 
One of them involves estimating the fraction of the events 
originating from the target cell walls and the other is an 
empirical overall normalization factor. 

To estimate the level of contamination from the target 
cell walls, events collected during the empty-target run 
period of the experiment are analyzed using the same 



process as those for the production run period. Only 
those events that fall within the target wall region for 
the empty target should be considered for the source of 
contamination. This is because even though there was 
no liquid hydrogen in the target, it was still filled with 
cold hydrogen gas. So, for this estimation only events 
within ±0.5 cm of the target wall region are selected. 
The correction is then calculated by taking the ratio of 
events within this target region from the empty target 
runs to those from the production run normalized to the 
total charge, p, collected during the run periods, 



R 



^*empty target /-^production 
production Pcmpty target 



(25) 



The average contamination is approximately 1% — 1.9% 
depending on the W kinematic bin. This ratio is then ap- 
plied as a correction factor to the measured cross section 
a = <r meas (l — R). Here, a is the corrected cross section 
and (J m eas is the measured cross section for a particular 
bin in W . 

The second correction (the empirical overall normaliza- 
tion factor) comes from comparing the measured ep — > ep 
elastic and the ep —> epir cross sections in the A (1232) 
resonance region (W — 1.23 GeV) to previously mea- 
sured values pH [28l EU |32] . The measured elastic scat- 
tering cross section from this experimental data were 
compared to the known cross section values [31) where 
both the electron and the proton were detected in the 
final state. A deviation of ~ 11% from the known cross 
section values is observed. 

This deviation of ~ 11% from the known elastic 
electron-proton scattering cross section includes the inef- 
ficiencies associated with the proton detection in CLAS 

mum. 

To account for this discrepancy, an overall normaliza- 
tion factor of Reiastic = 0.89 is applied to the ep — > epir 
differential cross section for every kinematic bin. An as- 
sociated systematic uncertainty of ±5% is applied. Af- 
ter this correction is applied, the measured ep —> epir 
cross sections for the A(1232) resonance region, W = 
1.23 ± 0.01 GeV, are in agreement with previous mea- 
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surements [2U [251 [32] to within 5% on average. Fig 
shows the result of this correction for a few kinematic 
bins in the A(1232) resonance region. 

Since the threshold region of interest for this experi- 
ment is sandwiched between the elastic and the A (1232) 
resonance region and the results in these two regions 
are consistent with previous measurements after applying 
this overall normalization factor, we believe this proce- 
dure is justified. This correction to the cross section also 
includes any detector inefficiencies and, as such, these 
inefficiencies will not be accounted for separately. 



VII. SYSTEMATIC STUDIES 

To determine the systematic uncertainties in the anal- 
ysis, the parameters of the likely sources of those uncer- 
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FIG. 12. The differential cross section ep — > epn for the 
A(1232) resonance region, W = 1.23 ± 0.01 GeV, for typical 
kinematic bins. The squares are the measured cross sections 
after applying the normalization correction factor (see text for 
details). The dashed curves are from Ref. and the dashed- 
dot curves are from the MAID2007 model. The corrected 
values agree with the two curves to within 5% on average. 



tainties are varied within reasonable bounds and the sen- 
sitivity of the final result is checked against this variation. 
A summary of the systematic uncertainties averaged over 
the kinematic bins of interest is shown in Table [H] 

The electron and proton identification cuts, the elec- 
tron fiducial cuts, the vertex cuts and the target cell cor- 
rection cuts provide small contributions to the overall 
systematic uncertainties. 

The electron EC sampling fraction cuts were varied 
from 3cr to 3.5a and the extracted structure functions 
changed by about 0.4% on average. The parameters for 
the electron fiducial cuts were similarly varied by about 
10% and the structure functions changed by about 1% 
on average. As such, a systematic uncertainty of 0.4% 
and 1% was assigned to these sources. 

The At cuts to select the protons were varied from 
3.5cr to 4cr and a variation of about 1.1% on average was 
observed on the extracted structure functions, which was 
assigned as the systematic uncertainty associated with 
this source. The variations in the fiducial cuts for the 
proton had a negligible effect on the structure functions. 

The vertex cuts were reduced by 5% and a variation 
of about 0.1% on average was observed on the extracted 
structure functions. So, a systematic uncertainty of 0.1% 
was assigned to this source. The structure functions are 
compared before and after applying the target cell cor- 
rections. A variation of about 1% is observed and this 
value was assigned as a source of systematic uncertainty. 

The major sources of systematic uncertainty are the 
Bethe-Heitler background subtraction, the missing mass 
squared cut to select the neutral pions, the elastic nor- 
malization corrections and the model dependence of the 
acceptance and radiative corrections. 



Source 


Estimate % 


e~ EC sampling fraction cuts 


0.4 


e~ fiducial cuts 


1 


p At cuts 


1.1 


Vertex cuts 


0.1 


Background subtraction cuts 


8 


7T° M x CUt 


3 


Target cell correction 


1 


Elastic normalization correction 


5 


Acceptance and radiative correction 


4 


Total 


10.8 



TABLE II. The average systematic uncertainties for the dif- 
ferential cross sections from various sources and the corre- 
sponding criteria. The final quoted systematic uncertainty, 
obtained by adding the different systematic uncertainties from 
each source in quadrature, is about 10.8%. 



There are residual Bethe-Heitler events that escape 
the elliptical Bethe-Heitler cuts. These events peak at 
M\ = 0, which have to be included in the overall fit. 
A Gaussian distribution was assumed for both the w° 
and the remaining Bethe-Heitler events. The pions are 
modeled by a Gaussian distribution near the expected 
pion mass and the Bethe-Heitler events are modeled by 



a Gaussian whose peak is at M\ = 0. This accounts for 
much of the tail in Figs. |8(b)| and |9| The resolution for 
M\ for the Bethe-Heitler and the pion distributions is 
expected to be similar because of the same kinematics 
of the detected electron and the proton. The Gaussian 
fit for the Bethe-Heitler is obtained, which is then sub- 
tracted to yield the pions. 

To see the effect of the background subtraction, the 
structure functions were compared with and without the 
application of the Bethe-Heitler background subtraction 
cuts. The structure functions changed by about 8% on 
average and this was used as a systematic uncertainty for 
this procedure. 

The missing mass squared cut was varied from 3<r to 
4ct and this resulted in a change of about 3% on average 
in the extracted structure functions. 

The systematic uncertainty on the elastic normaliza- 
tion correction of ±5% was obtained by looking at the 
difference between the extracted structure functions be- 
fore and after applying the correction factor to the data. 
The structure functions varied by about 5% on average. 

Additionally, a ±4% systematic uncertainty is assigned 
on the model dependence of the acceptance and radiative 
corrections based on previous analyses [TH1 [Ml HI] . 

The total average systematic uncertainty, obtained 
by adding the individual contributions in quadrature is 
10.8%. 
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VIII. DIFFERENTIAL CROSS SECTIONS AND 
STRUCTURE FUNCTIONS 



at other kinematics is similar to the MAID2007 model. 



The kinematic coverage of the experiment spans over 
W from 1.08 to 1.16 GeV and Q 2 from 2 to 4.5 GeV 2 . 
The reduced differential cross section for the reaction is 
computed for each kinematic bin. The cross sections are 
reported at the center of each kinematic bin. Fig. |l3| 
shows the differential cross section for some of the kine- 
matic bins near threshold as a function of (f>*. The pre- 
dictions from LCSR Q3J, MAID2007 [2E] and SAID [33] 
are shown for comparison. 



Using Eq. (161, the differential cross section is fitted to 
extract the structure functions <jt + ^l-, &tt and ct.lt- 
The result of the fit is shown as the solid curve in Fig. fl3] 
The reduced \ 2 f° r the fit is calculated using \ 2 = x\jv, 
where v is the number of degrees of freedom calculated 
for each W, Q 2 , and cos #* bin (i.e., v = 6 data points —3 
fit parameters = 3), and \o is the unnormalized goodness 
of fit. The averaged x 2 of the fits is 0.9. 

The extracted structure functions <jt + ecl , &tt and 
olt are shown in Figs. [Wj [15] and [l6j respectively, as 
a function of cos0* for W = 1.08 - 1.16 GeV and 
Q 2 = 2.0—4.5 GeV 2 . The data points are shown with sta- 
tistical error bars only and the size of the systematic er- 
rors is shown as the gray boxes. Predictions from LCSR, 
MAID2007, and SAID are also included for a T + ea L and 
<jlt- Since the LCSR does not include any utt contri- 
butions in the calculations, they are not shown. 

The structure function ctt + £Ol (Fig. 14) is generally 
in good agreement with the MAID2007 predictions but 
there is some discrepancy for W = 1.09 GeV at high 
cos 9*. This discrepancy is reduced for higher W bins. 
The results disagree with the LCSR predictions, espe- 
cially for those bins away from threshold (W > 1.09 
GeV). This disagreement is also apparent for low Q 2 
bins. As one moves closer to threshold and at high Q 2 , 
the agreement is quite good, especially at backward an- 
gles cos 6** — > — 1. The LCSRs have been calculated and 
tuned especially for the threshold region at high Q 2 and 
thus, there exists a strong disagreement at higher W and 
low Q 2 bins. The predictions from SAID strongly dis- 
agree for the first W bin and low Q 2 bins, but converge 
toward the MAID2007 predictions for higher W and Q 2 . 



The structure function utt (Fig. 15) results are in 



good agreement with the SAID and MAID2007 predic- 
tions for low W and high Q 2 but disagree at high W and 
low Q 2 bins. Most of the values are close to zero for all 
W. The LCSR predictions assume only s-wave contri- 
butions to the cross section from this structure function. 
The <i-wave contribution to the total cross sections in 
SAID range from to 0.001 /ib for the near threshold 
bins m. 



The structure function ctlt (Fig. 16) also shows good 
agreement with the MAID2007 and LCSR predictions 
for high Q 2 and low W, but there is some discrepancy at 
other kinematics. The SAID prediction has a large dis- 
agreement at low W and Q 2 , but the level of agreement 



IX. S-WAVE MULTIPOLES AND 
GENERALIZED FORM FACTORS 

In order to compare with the calculated generalized 
form factors of Ref. [14], one must extract the s-wave 
multipole amplitudes from the measured cross sections. 
First, the structure functions are written in terms of the 
helicity amplitudes Hi . The helicity amplitudes are func- 
tions defined by transitions between eigenstates of the 
helicities of the nucleon and the virtual photon [TB] . The 
helicity amplitudes are then expanded in terms of the 
multipole amplitudes. 

The structure functions are related to the helicity am- 
plitudes iii, 2 ,... 6 (VF,Q 2 ,cos(9;) by: 



VT = lm\ 



\H 2 \ 2 + \H 3 \ 2 + \H4 2 ), 



a L = \H 5 \ 2 + \H 6 \ 2 , 

1 



(26) 

(27) 
(28) 



V2 



ReHHi - Hi)H* + (H 2 + H 3 )H* 6 ]. (29) 



The analysis of the data is based on the following ex- 
pansion of the helicity amplitudes over multipole ampli- 
tudes (see, for example, [35]): 



Hi = sin 0: cos 9 f ]T(i? ;+ - 
[P('(co S 6* w )-P(' +1 (cosd;)], 

[iftcoBO-^+i(ca<)]. 



1 

[P/'tcos^ + P/i^cos^)] 



H 3 = — sine; sin f ^(B i+ + B (J+1) _) 



H 4 - V2sin f ]T(A i+ + A (l+1) _) 
[P/(cos^) + P/ +1 (cosC)], 



cos^£(Z + l)(S ;+ + S (m) _) 



Q 

— r cos „ 

|q*l 2 

[P/(cos^)-P/ +1 (cos^)], 

Q ■ 

- — r sin 
|q*l 2 



(30) 
(31) 

(32) 
(33) 
(34) 



11 ^sin^^(Z + l)(S i+ -S (i+1) _) 

[P/(cos^) + P/ +1 (cos^)]. (35) 



Here, P[ l+1 (cos8^) and P[' l+1 (cos8^) are the first and 
second derivatives of the Legendre polynomials, respec- 
tively, and q* is the virtual photon 3-momentum in the 
cm. system. Also, 



a,.+ = \[(i + 2)£j+ + m+] , 



B, 



1+ 



E, 



1+ 



Mi 



1+' 



(36) 
(37) 



13 





FIG. 13. The differential cross sections in ^tb/sr for a few kinematic bins near threshold as a function of (f>%. Experimental 
points (squares) are shown with statistical uncertainties only. The size of the estimated systematic uncertainties is shown in gray 
boxes below. The predictions from LCSR, MAID2007 and SAID are shown as dashed, dashed-dotted and dashed-double-dotted 
curves, respectively. The fit to the distributions is shown as a solid curve. See Sec. |VHll for details. 



-4, 



1 



i(i+D- = 2 [V + 2 ) M ('+D- - lE (l+i)-] > (38) 
B (i+i)- = E (i+i)- + M (i+\)-- (39) 

The strong cos 6* -dependence of the structure function 
<jt + £&l, and the nonzero values of <jlt found in the ex- 



periment (see Figs. 14 and 16 ) show that higher multipolc 
amplitudes should be taken into account in addition to 
the s-wave amplitudes Eq + and 5*0+ at all W . Our un- 
derstanding of the high-wave multipoles, which should 
be included in this analysis, was based on the results of 
the analysis of CLAS data [Ml US] performed in Ref. [32] 
using the unitary isobar model (UIM) and dispersion re- 
lations (DR). These data are on the j*p -+ 7r + n [25] and 
7*p — » 7r°p [M] cross sections in a similar range of Q 2 but 
in a significantly wider energy range, which start from 
W = 1.15 and 1.11 GeV, respectively. The precision in 
the present experimental results near threshold is much 
better than the precision in Refs. [Ml US]. However, the 
results of their analysis are useful to study the p- and d- 



wave contributions, which are determined mainly by the 
A(1232)P 33 , iV(1440)P u , and N(1520)D 13 resonances. 

According to the results of the analysis [32] at W = 
1.09 to 1.15 GeV, there are large p-wave contributions 
related to the A(1232)P 33 and iV(1440)Pii. The d-wave 
contributions are negligibly small for the following rea- 
sons: (i) near threshold, the c?-wave multipolc amplitudes 
are suppressed compared to the p-wave amplitudes by 
the additional kincmatical factor p*; (ii) at the values 
of Q 2 investigated in this experiment, the contribution 
of the A^(1520)£>i 3 to the corresponding multipole am- 
plitudes is significantly smaller than the contributions 
of the A(1232)P 33 and iV(1440)Pn to the p-wave mul- 
tipole amplitudes; (hi) in contrast with the A(1232)P 33 
and AT(1440)P U , the width of the Ar(1520)D 13 is signif- 
icantly smaller than the difference between the mass of 
the resonance and total energy at the threshold. There- 
fore, in our analysis only multipole amplitudes P + 1 $0+ 1 
Mi±, S\±, and E\+ were included. 
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FIG. 14. The structure function or + £(Tl as a function of 
cos6* in jitb/sr for W = 1.08 - 1.16 GeV and Q 2 = 2.0 - 4.5 
GeV 2 . Predictions from LCSR that include only s-wave 
contribution (dashed), MAID2007 (dashed-dot), and SAID 
(dashed-double-dot) are shown. The error bars represent sta- 
tistical uncertainties only and the estimated systematic uncer- 
tainties are shown as gray boxes. The solid curve corresponds 
to t he r esults obtained from the fit to the cross sections (see 
Sec. |lX| for details). The values of Q 2 (on top of the panels) 
and W (on the right side of the panels) are the central values 
of the bins. 



FIG. 15. The structure function gtt as a function of cos#£ 
in ^b/sr for W = 1.08 - 1.16 GeV and Q 2 = 2.0 - 4.5 GeV 2 . 
Predictions from MAID2007 (dashed-dot) and SAID (dashed- 
double-dot) are shown. The LCSR predictions do not include 
any utt contributions, so they are not shown. The error 
bars represent statistical uncertainties only and the estimated 
systematic uncertainties are shown as gray boxes. The solid 
curve corresponds to t he r esults obtained from the fit to the 

for details). The values of Q 2 (on 



IX 



cross sections (see Sec. 
top of the panels) and W (on the right side of the panels) 
are the central values of the bins. The horizontal line at zero 
serves as a visual aid. 



The data were fitted simultaneously at W = 1.09, 1.11, 
1.13 and 1.15 GeV with statistical and systematic uncer- 
tainties added in quadrature for each point. The am- 
plitudes were parametrized according to their threshold 
behavior and the results of the analysis in Ref. [52] . 

Due to the Watson theorem the imaginary parts 
of the multipole amplitudes below the 2n production 
threshold are related to their real parts as ImM. = 
ReMta,n(5^ N ), where M denotes E[±, M( ± or S[ ± am- 
plitudes, and I is the total isotopic spin of the irN sys- 
tem. Near threshold S^ N ~ p* 2i+1 , and the imaginary 
parts of the multipole amplitudes are suppressed com- 
pared to their real parts. Therefore, in the analysis, only 



the real parts of the amplitudes were kept. These ampli- 
tudes were parameterized as follows: Eq + ,Sq+ ~ const, 
M 1± , Si±, and E 1+ ~ p*. 

In the fitting procedure, the amplitudes Eq + , So+ and 
Mi± were fitted without any restrictions. The relatively 
small amplitudes S±± and E\+ were fitted within ranges 
found from the results of the analysis of the data [2H 
using the UIM and DR in Ref. 32J. It should be men- 
tioned that the results for the M\± contributions ob- 
tained in our fit of the 7*p — > ir°p cross sections near 
threshold are consistent with those of Ref. [35] obtained 
in the analysis of significantly larger range over W. The 
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FIG. 16. The structure function glt as a function of cos#^ 
in ^b/sr for W = 1.08 - 1.16 GeV and Q 2 = 2.0 - 4.5 GeV 2 . 
Predictions from LCSR that include only s-wave contribu- 
tion (dashed), MAID2007 (dashed-dot), and SAID (dashed- 
double-dot) are shown. The error bars represent statistical 
uncertainties only and the estimated systematic uncertainties 
are shown as gray boxes. The solid curve corresponds to the 
results obtained from the fit to the cross sections (see Sec. IX 
for details). The values of Q 2 (on top of the panels) and W 
(on the right side of the panels) are the central values of the 
bins. The horizontal line at zero serves as a visual aid. 



overall average \ 2 P er degree of freedom for the fit is 
approximately one. 

The obtained results for the structure functions are 
plotted in Figs. 14p6 as solid curves. It can be seen that 
the multipole amplitudes Eq + , <So+, M±±, Si±, and E\+ 
parametrized in the way discussed above represent the 
data very well at all W . The obtained results for E$ + and 
5*0+ are presented in Fig. 17 These multipolcs have been 



normalized to the dipole formula Gd{Q ) = ^1 + ~ n 

Fig. [18] shows the extracted generalized form factors, 
Gi and G2, as a function of Q 2 . The error bars on 
the points include statistical and systematic uncertain- 
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FIG. 17. The s-wave multipoles (a 



and (b) So+ normal- 
ized to the dipole formula Go are plotted as a function of Q 2 . 
The error bars include statistical and systematic uncertainties 
added in quadrature. The size of the estimated systematic un- 
certainties are shown in the bottom. The LCSR based model 
predictions and the LET predictions are also shown as curves. 
The horizontal line at zero serves as a visual aid. 



ties added in quadrature. The size of the estimated sys- 
tematic uncertainties is shown separately at the bottom 
of the plots, which assumes all systematic errors for all 
the data points to be entirely uncorrelated (10.8%). The 
LET 7] predictions are shown as dash-dotted curves. 

The plots also show LCSR predictions [H] as solid and 
dashed curves. Braun et al. have tried to minimize the 
uncertainties in their LCSR based model calculations by 
including electromagnetic form factor values known from 
experiment. These calculations are shown as solid curves 
in the figure. The "pure" LCSR based models are cal- 
culations where all the form factors are obtained entirely 
from theoretical calculations and the uncertainties have 
not been minimized. These are shown as dashed curves 
in the figure. The difference between these two curves 
can essentially be treated as the overall uncertainty in 
their predictions. 
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FIG. 18. The generalized form factors (a) G\ p and 
normalized to the dipole formula Gd are plotted as a func 
tion of Q 2 . The error bars include statistical and systematic 
uncertainties added in quadrature. The size of the estimated 
systematic uncertainties are shown in the bottom. The LCSR 
based model predictions and the LET predictions are also 
shown as curves. The horizontal line at zero serves as a vi- 
sual aid. 



X. DISCUSSION 

The results for the E + multipole and G\ p are in good 
agreement with the LCSR predictions. The extracted 
E 0+ values deviate significantly from the LET predic- 
tions over the entire Q 2 range even though the extracted 

G\ p values are not too far off from the LET predictions. 
This is because the LET calculations for Eq + only de- 
pend on G\ p (Eq. whereas the LCSR calculations 
include contributions from both G\ p and G\ p (Eq. &2i). 

' — "o 

The overall trends of increasing Eq+ and decreasing G\ p 
are similar to these two predictions, but the deviation of 

the extracted values for G\ p from the LET predictions 
becomes much more apparent at Q 2 > 3 GeV 2 . 

One can observe a discrepancy of our results for the 



So+ multipole and G\ p from the LCSR predictions. The 
results are closer to the LET predictions but are not en- 
tirely consistent for all Q 2 . 

The uncertainty in the LCSR predictions for the Sq+ 



multipole and G\ p is much bigger than for Eq + and 

G\ p . In the chiral limit approximation, m, v — > 0, the 
Pauli form factor F2(Q 2 ), which is the primary contribu- 
tor to the calculations of So+ and G\ p , is not reproduced 
very well. Also, the LCSR calculations exist in lead- 
ing order only and do not include next-to-leading order 
(NLO) corrections. The NLO corrections are expected to 
be large. Additionally, the LCSR predictions contain ap- 
proximations and were not expected to have an accuracy 
of better than 20% [13] . 

Furthermore, the LCSR predictions do not include ef- 
fects from terms proportional to the pion mass. In the 
Q 2 region of this experiment, the predictions indicate a 
suppression of the Sq+ multipole [HJ and this multipole 
is very sensitive to corrections of all kinds, including the 
pion mass corrections. In the LET predictions, some pion 
mass corrections have been included [7 . This may also 
explain the discrepancy between the predictions and the 
extracted results for Sq+ and G\ p . 

Due to these theoretical uncertainties, the predictions 

of the magnitude of Sq+ and G\ p /Gd, and where they 
cross zero, differs for the two methods of calculation. The 
experimental results indicate that this sign change for 

G\ p /G D occurs at Q 2 > 4 GeV 2 rather than at the 
LCSR prediction of around 2.2 GeV 2 or 3.5 GeV 2 . 

The results of the structure functions, Figs. 14fT6 indi- 
cate a significant contribution of the p-wave in the near 
threshold region as indicated by the almost linear de- 
pendence of the ot + £<Jl as a function of cos0*. This 
contribution increases as one moves away from threshold 
to higher W (e.g., see Fig. 16). This is highly under- 



estimated in the overall LCSR predictions for the struc- 
ture functions and cross section calculations. Their pre- 
dictions are tuned to include mostly s-wave and very little 
p-wave contribution very close to threshold at high Q 2 . 
This also explains the good agreement of the extracted 
Eq+ and G\ p to their predictions but the strong dis- 
agreement of the So+, G\ V i the cross sections and the 
structure functions. 

The extracted generalized form factors, G\ p and 

p , show a faster fall off than the dipole form. This 
suggests a broadening of the spatial distribution of the 
correlated pion-nucleon system. It suggests that the cor- 
related pion-nucleon system is broader than the bare nu- 
cleon itself because the bare nucleon follows the dipole 
form factor. 

The results for G\ p show similar trends to the previ- 
ously extracted G\ +n [TS]. In comparison, the former is 
about 30% higher in magnitude while the overall behav- 
ior as a function of Q 2 is similar. There are no results for 
Gl 71 for comparison. However, the generalized form fac- 
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tor results for the Tr a p channel provide strong constraints 
on chiral aspects of the nucleon structure and the validity 
of the LETs at high Q 2 . 
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